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Abstract 

We present heterodyne detected transient grating measurements on water filled Vycor 7930 in 
the range of temperatm'e 20 — 90 °C. This experimental investigation enables to measure the 
acoustic propagation, the average density variation due the liquid flow and the thermal diffusion in 
this water filled nano-porous material. The data have been analyzed with the model of Pecker and 
Deresiewicz which is an extension of Biot model to account for the thermal effects. In the whole 
temperature range the data are qualitatively described by this hydrodynamic model that enables 
a meaningful insight of the different dynamic phenomena. The data analysis proves that the signal 
in the intermediate and long time-scale can be mainly addressed to the water dynamics inside the 
pores. We proved the existence of a peculiar interplay between the mass and the heat transport 
that produces a flow and back-flow process inside the nano-pores. During this process the solid 
and liquid dynamics have opposite phase as predicted by the Biot theory for the slow diffusive 
wave. Nevertheless, our experimental results confirm that transport of elastic energy (i.e. acoustic 
propagation), heat (i.e. thermal diffusion) and mass (i.e. liquid flow) in a liquid filled porous glass 
can be described according to hydrodynamic laws in spite of nanometric dimension of the pores. 
The data fitting, based on the hydrodynamic model, enables the extraction of several parameters 
of the water- Vycor system, even if some discrepancies appear when they are compared with values 
reported in the literature. 

PACS numbers: 78.47.jj, 47.61.-k, 62.80.+f, 47.56. +r 
Keywords: 



I. INTRODUCTION 

The study of transport phenomena in heterogenous media is a fundamental issue of the 
material science l|, |2|, its relevance spans from the basic physics (i.e. the proper definition 
of the equations of motion describing the transport processes) to the more recent techno- 
logical applications. During the recent years, these studies have to face the scaling down of 
heterogeneity towards the nano-metric dimension. 

Between the numerous random heterogenous media the two phase solid-liquid materials 
represent one of the most investigated. A large impulse, to these studies, has been given by 
the petroleum industry aimed at understanding the transport phenomena in sediment sands 
and porous rocks. In the basic research, probably the most studied materials are the solid 
porous matrix filled by a molecular liquid (e.g. liquid filled porous glasses), since typically 
they are well parameterized systems. 

A part from the electric processes, three main transport phenomena are generally to 
be considered in the solid-liquid heterogenous materials: propagation of acoustic waves, 
thermal diffusion and viscous flow (i.e. transport of elastic energy, heat and mass). These 
phenomena are generally characterized by very different time/frequency scales so that often 
they have been considered as quasi-independent processes, both from the experimental and 
theoretical point of view. Moreover, they can be surprisingly described by relatively simple 
hydrodynamic models. 

The propagation of sound in the liquid-filled porous materials has been described by phe- 
nomenological models based on hydrodynamic/elastic equations for the liquid/solid phases. 
The first model has been introduced by M. A. Biot in an important theoretical work in 
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J]. An extraordinary result of the Biot theory is the prediction of a new slow longi- 



tudinal wave besides the usual longitudinal and transverse waves. This wave has a velocity 



lower than the one of the bulk liquid 



B 



6|. In a recent experimental work the validity of this 



model in the high frequency range has been investigated |7|, l8[ , measuring the hypersonic 
sound propagation in a nano-porous media. 

The viscous flow of liquids through porous media has always been a subject of intense 
study |9|. Many studies have been concerned with the check of the validity of the macroscopic 
law of Darcy. Although Darcy's law was determined phenomenologically, it is a direct 
result of the Navier-Stokes hydrodynamic equations 10||, thus checking its validity is a 



straightforward analysis of the hydro dynamical character of liquid flowing. Darcy's law is 



shown to be valid for a wide range of porous media from porous rocks, sands 



3 and glass 



beads [llj whose heterogeneities are micrometric, up to porous glasses like Vycor which have 
heterogeneities of the order of few nanometers 12|-|l6l|. Recently, many experimental and 



simulation works aimed at understanding to what extent the Navier-Stokes equation can 



describe the liquid flow at decreasing of the confinement size [ij, Il7l-ll9|. 



The thermal conductivity of heterogenous media is a complex problem which has to be 
characterized in a general theoretical framework [ij. The prediction of the effective per- 
meability of saturated porous media remains, despite of many experimental and theoretical 
works, an unsolved problem in heat transfer science. In particular this quantity depends, 
as well as on macroscopic parameters of the two constituent phases, also on the particular 
morphology of the porous material which is in general experimentally not accessible. Nev- 
ertheless, the theory is able to fix the limits of the effective conductivity and give several 
approaches for calculating this quantity in particular kind of porous media |20| . 

Though the transport phenomena in liquid-filled porous glasses has been previously stud- 
ied in the literature, many basic questions remain open. In our opinion, one of the more 
relevant is at which extent the hydrodynamic models are valid, especially when the media 
heterogeneity scale down versus molecular length scales. 

We think that other experimental researches could infer new information, especially with 
techniques poorly applied in this field, like the time resolved spectroscopy experiments. 
Transient grating (TG) experiments 2l|, |22[ are powerful tools for investigate the relaxation 



dy„a„,. of complex liquid. ^^1. but only few p.evio.. ™.ental wo.k. util^ed the=e 



techniques to investigate the porous glass samples |7|, |8|, |26|, |27| . A particular kind of TG 
experiment offers the possibility to study, at one time, the acoustic propagation, the liquid 
flow and the thermal diffusion processes. The very broad time window covered by this 
experiment, typically from 10^^ to 10^^ s, gives access to a dynamic range hardly explored 
by other methods. 

In this paper we present heterodyne detected transient grating (HD-TG) measurements 
on water filled Vycor nano-porous glass. The measured signal shows aspects related to the 
acoustic waves propagation, to the viscous liquid flow through pores, and, lastly, to the 
thermal diffusion. The transient grating experiment reported here is shown to be a powerful 
tool to measure the transport processes in liquid filled porous glasses. 



To a naly ze the data we employ an hydro dynamic model introduced by Pecker and Dere- 
siewicz 28| which is a extention of the Biot's theory |3|, y] to include the thermal effects. 



The model enables a safe and meaningful addressing of the different contributions present 
in the HD-TG signal, describing them on the base of few parameters characterizing the 
heterogenous system. 

The present work is divided as follows. In Sec. [TTl we report an overview of the TG 
experimental technique describing the involved excitation mechanisms and the measured 
dynamics with probing process. In Sec. IIIII we present the experimental aspects of the 
work, laser systems, experiment optical set-up and the sample preparation. Section IIVI is 
devoted to the presentation of the data while the subsequent one to the introduction of the 
theoretical model used to analyze these data. Finally, in the last section we show the data 
analysis and we discuss the obtained results. 

II. TRANSIENT GRATING EXPERIMENTS 

In a TG experiment, two infrared laser pulses, obtained dividing a single pulsed laser 
beam, interfere within the sample and produce an impulsive spatially periodic variation of 
the dielectric constant. The spatial modulation is characterized by a wave vector q which 
is given by the difference of the two pump wave vectors ki — k2 (see Fig. [1]). Its modulus 
is g = 47rsin(6'ex)/Aea;, where \ex and 9ex are the wavelength and the incidence angle of the 
exciting pumps, respectively. The relaxation toward equilibrium of the induced modulation 
can be probed by measuring the Bragg scattered intensity of a second cw laser beam. The 
time evolution of the diffracted signal supplies information about the dynamic of the relaxing 
TG and, consequently, on the dynamical properties of the studied sample. 



TG experiments fall within the framework of the four-wave-mixing theory 2l|, |25|, l29|, l30 |. 
Under a few approximations, in particular, if the laser pulses do not have any electronic 
resonance with the material, it can be proved that the TG experiment can be divided in two 
separated processes: excitation and probing. Moreover, when the heterodyne detection is 
employed, the signal turns out to be directly proportional to the dielectric constant change 



induced by the pumps, 5eij, or to the response function of the system, Rijki 25 1 



S (g, t) ex 6eij{q, t) ex Rijki{q, t) (1) 




FIG. 1: Schematic drawing of a transient grating experiment. Two excitation pulses, Ef^ and Ef^ 
induce an impulsive spatial modulation of the dielectric constant with step A. The relaxation of 
the induced transient modulation is probed by the Bragg scattering of a third beam, Ep. 

where the cartesian indexes i and j define the polarizations of the diffracted and the probe 
fields, k and / define the polarizations of the pump fields. Both 6eij and Rijki are spatial 
Fourier components corresponding to the exchanged wave vector q, i. e. the grating wave 
vector. Hence, the heterodyne detected TG signal measures directly and linearly the re- 
laxation processes defined by the tensor components of the response function. By selecting 



different polarizations of the fields, different elements o 
probed. Generally, these elements can be different 
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the response function tensor are 



34| . In materials, where the coupling 



between translational and rotational degrees of freedom is strong, the birefringence effects 
can be relevant and the TG signal will depend on the field polarization. When, instead, all 
the birefringence contributions are negligible, the signal turn out to be independent on the 
field polarization and the induced dielectric constant tensor is simple I6e, where / is the 
identity tensor [25 1. 

Without birefringence effects and in bulk homogeneous materials, the transient grating 
is mainly induced by two effects. The electrostriction and the heating. The first effect 
arises from the interaction between the molecular dipoles induced by the pump fields and 
the same pump fields. The electrostriction produces a pressure grating and consequently 
a density grating launching two counter-propagating acoustic waves whose superposition 
makes a standing acoustic wave. The second exciting process arises from a weak absorption 
of the pump infrared radiation which is resonant with some vibrational states. The increasing 
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of energy due to the pump absorption, generally, quickly thermalizes by fast non radiative 
channels and builds up a temperature grating. This one produces a pressure grating and 
then the pressure grating produces a density grating via thermal expansion. Now, besides 
acoustic waves, we have more a constant density grating supported by the temperature 
grating which relaxes by thermal diffusion. As we shall see in section IIV[ in a heterogeneous 
material the excitation sources are the same but the effects can be more complex. 

III. EXPERIMENTAL PROCEDURES 

A detailed description of the lasers and the optical set-up of our TG experiment is reported 



in references (25|. Here we just recall the main aspects. A sketch of the TG optical set-up 
is reported in Fig. [2l The infrared pump pulses have a 1064 nm wavelength, a temporal 
length of 20 ps and repetition rate of 10 Hz. They are produced by an amplified regenerated 
oscillator (Nd-YAG EKSPLA PL2143). The typically used pump pulse energy was 5 mJ. 
The probing beam, instead, is a continuous-wave laser at 532 nm produced by a diode- 
pumped intracavity-doubled Nd-YVO (Verdi-Coherent). The two laser beams are collinearly 
sent to a phase grating (PG) to get the two pump pulses and the probe and local field beams. 
These are obtained taking the -|-1 and —1 diffraction orders of the infrared and green lasers. 
A couple of confocal achromatic lenses, ALl and AL2, collects all the four beams and focuses 
them on the sample. The phase grating directly supplies a probe at the right Bragg angle and 
a local field exactly coUinear with the scattered field and phase locked with the probe. The 
HD-TG signal is optically filtered and measured by a fast avalanche silicon photodiode with 
a bandwidth of 1 GHz (APD, Hamamatsu). The signal is then amplified by a DC-800 MHz 
AVTECH amplifier and recorded by a digital oscilloscope with a 7 GHz bandwidth and a 
20 Gs/s sampling rate (Tektronix). The instrumental function of our setup has a temporal 
full width half maximum of 1 ns and is mainly determined by the bandwidth of the detector 



and its amplifier 25 1. 



Vycor 7930 is an open cell porous glass which is produced by a spinodal demixing in a 



borate-rich and borate-poor glass and subsequent bleaching of the borate- rich phase (35|. It 
has nominal values of porosity and mean pore size of 28% and 4 nm respectively. The solid 
constituent of Vycor 7930 is the glass Vycor 7913 which is composed by 96% of silica, by 
3% of boron oxide and the remaining part mostly by aluminum oxide and zirconium oxide. 
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FIG. 2: Optical set-up and laser system for the TG experiment with optical heterodyne detection. 
M are the mirrors; CL^^ the cylindrical lenses; DM a dichroic mirror; PG is the phase grating; 
AL^^ the achromatic lenses and PD is the photodiode. 

Our sample has been purchased from "Advanced Glass and Ceramics" in a cylindrical form, 
with a diameter of 15 mm and 10 m,m, thickness. In order to remove all the absorbed 
organic elements we washed the sample in 30% water-hydrogen peroxide solution and then 
we heated it in a muffle up to 800 °C at a rate of 0.3 K/min. The sample has been kept at 
this temperature for several hours. Water has been obtained from a double-distilled water 
vial prepared for pharmaceutical purposes. The filled matrices were, then, placed in a thin 
teflon tube and closed between two circular quartz windows. The whole was inserted in a 
cylindrical copper cell (similar to the one reported in J36|). The holder was then stabilized 
in temperature with a stability of ±0.1 K. 



IV. EXPERIMENTAL RESULTS 



The HD-TG signals on water-fllled Vycor have been collect in the range of temperature 
20 — 90 °C at different wave vectors. The data did not show a dependence on the polariza- 
tions of beams and consequently birefringence contributions were negligible 2J] . Hence, we 
recorded the data in a single conflguration of beams polarization in which all the beams had 
a polarization perpendicular to the scattering plane as sketched in Fig. [H The data at three 
different temperatures, 90, 40 and 20 °C, are shown in Fig. [3] (following the usual TG signal 
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FIG. 3: HD-TG data on Vycor-water system at g = 1 fim~^ at three different temperatures. At 
sliort times, tlie data show damped acoustic oscillations, at middle times, a density rearrangement 
related to liquid ffow and finally, at long times, the data show a decay due to the thermal diffusion. 
The acoustic waves are practically induced only by the electrostrictive effect, while the remaining 
part of the signal by the thermal excitation. 



convention, tlie sign of tlie data is chosen positive for a negative change of density 25]). 

Before explaining the data features, some considerations have to be discussed. As just 
mentioned in section [Til in a TG experiment we have two exciting sources: the electrostrictive 
effect and the heating process. In the water filled Vycor, the first excitation is surely present 
in both materials, even if the electrostrictive strength could be different in the two materials, 
definitely we expect that this excitation source launches the acoustic waves. Differently, 
the heating should be much weaker in the solid part than in the liquid water or at the 
pore surfaces. In fact, the liquid water has a no negligible absorption coefficient at the 
pump wavelength and also the pore surfaces present a weak infrared absorption due to the 
presence of the silanol and boranol groups (Si-OH and B-OH). Nevertheless, on a very fast 
time scale, in the illuminated areas, all the sample components (solid part of the matrix, 
pore surfaces and liquid water) reach the same temperature. This makes applicable the 
assumption of local thermal equilibrium, i.e. we are supposing that the liquid and matrix 
temperatures are always locally equilibrated. Kaviany ([33] pp. 120-121) reports criteria 
of validity for the local thermal equilibrium approximation which are well fulfilled for the 



water- Vycor system. Namely, the approximation of local thermal equilibrium requires that 
the time scale associated with the inter-pore heat transfer must be much smaller than the 
time scale associated with the macroscopic heat transfer experimentally measured. The heat 
transfer between two pores occurs on a nanometric length scale while the macroscopic length 
scale is given by the spacing of the induced grating, A = 27r/g = 6.28 fim, for q = 1 fim~^. 
Moreover, the interphase heat transfer is a very efficient process thanks the very large specific 
pore surface area. We expect this process to evolve on shorter time scales than the time 
scales measured in the HD-TG signal. So, also in our heterogenous sample the pump heating 
produces an uniform temperature grating as in a bulk homogeneous sample. 

Now, other important aspects have to be considered to understand the HD-TG signals: 
1) the thermal expansion coefficient of the matrix is very low compared to the water one; 2) 
the bulk modulus of the matrix is higher than the water one, the matrix is stiff. Since the 
thermal expansion of water would be greater than that of the matrix, at first the expansion 
of the liquid is reduced by the stiffness of the matrix which exerts a pressure on the liquid. At 
later times, the liquid flows through the pores to nullify the matrix pressure which is at flrst 
positive and then negative, as it will be proved later using a data simulation based on the 
hydrodynamic model. So, the liquid can practically expand only via the flowing. This implies 
also that the acoustic oscillations induced by the temperature grating are characterized by 
a very low amplitude. 

For the aforementioned reasons, the data show, at short times, damped acoustic oscil- 
lations induced almost only by the electrostrictive effect, at intermediate times, a density 
rearrangement related to liquid flow inside the pores. 

From preliminary flts of the intermediate and long parts of the signals, we obtain that 
the rise is well reproduced by a single exponential whose time constant, Ty, shows a q~^ 
dependence. It is, therefore, a diffusive mode and it turns to be the Biot slow wave in diffusive 
regime, as it will shown later using the hydrodynamic model. The time constant, then, is 
related to the constant diffusion of the Biot slow wave (hydraulic diffusivity coefficient) 
which is connected to Vycor permeability and liquid viscosity. Clearly, in this case the 
temperature grating plays an important role, it is the pumping source for the liquid motion. 
In Fig. |4]we report the inverse of constant times of the exponential signal rise, Ty = 1/Ty, 
and the exponential signal fall, Tt = ^/tt, as a function of the g^ for the temperature fixed 
at 30 °C. The probed wave vectors were q = 0.63, 1.00, 1.39, 1.76, 2.15, and 2.51 ^m~^. 
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FIG. 4: Decay constants of the two exponentials describing the long part of the signal as a function 
of q^ at the fixed temperature of 30 °C. Vy = 1/ry and Fy = l/rr are, respectively, the decay 
constants of the exponential rise and fall of the signal. The linear behavior confirms the diffusive 
character of the viscous flow mode and of the thermal diffusion. 

The linear behavior of the time constants confirms the diffusive character of the modes. 
V. THEORETICAL BACKGROUND 



An exphcit function of the TG signal can be obtained developing the dielectric function be 
through the hydrodynamic variables which, in the absence of birefringence effects, are only 
the density and the temperature. According to a first-order approximation, the dielectric 



constant change in a homogeneous bulk material is J 



5T(t), 



(2) 



Now we have to consider that our sample is composed by two interconnected materials, 
Vycor and water. The dielectric constant change will depend on two densities and two 
temperatures: 



be{t) 



de 
dpi 



Tm,Tl,p2 \ ^-^ / Tm,Tl,pi 



m{t) 



\9^Jt„p„pJ'^"'^^^^\^^^ 



(3) 



Tm,pl,p2 

where pi = (1 —(p)ps and p2 = (pPi are the average densities of matrix and water respectively, 
in which (p is the porosity of Vycor, p^ the density of solid and pi the water density. T^ and T/ 
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are the matrix and water temperatures respectively. The dielectric constant change induced 
directly by the temperature variation, in Eq. (|2]), is generally neglected in the TG signal 
analysis because for most of liquids {de/dT)p6T is much lower than {de/dp)TSp. Water is 
an example for which this condition is not satisfied, and the TG signal depends also on this 
contribution 2J]. Thus, we can neglect the matrix temperature contribution, but not the 



water temperature contribution that could have an important weight: 

de 



6e{t) 



dpi 



Spiit) 



|^)wt) + (^KTKt) 



The previous equation can be rewritten in the following way: 



^^'"^'H 



Ai25pi{t) + 5p2{t) + E^dTiit) 



where we have defined 



A 



12 



{de/dpi 



and E 



Toide/dTi) 



(4) 



(5) 



(6) 



{de/dp2) p2 {de/dp2) 

Ayi and E are amplitude parameters which define the relative weight between the three 
different contributions in the final signal. E is defined similarly in 2J], but its value could be 
different from that measured for the bulk water because of the strong interactions between 
water molecules and the hydrophilic surfaces of Vycor pores. This interaction strongly 
modifies the hydrogen bonds among water molecules and it could change the value of E. 
Thus, both the parameters, A\2 and E must be valued from the data fitting. 



A. Hydro dynamic model 



The time evolution of the two densities and the two temperatures can be obtained from 
the thermo-poroelastic model introduced in 1973 by Pecker and Deresiewicz (P-D) 



This is an extension of the well known Biot model on the wave propagation on poroelastic 
system to account for the temperature effects. An earlier attempt to include the temperature 
was made by Zolotarev in 1965 39|, but here strong approximations were introduced. The 
model of Pecker and Deresiewicz has been later deeply revisited as linearized approximation 
of more complex models by Gajo 40|] and Youssef 4l|. We shall follow here the Pecker and 
Deresiewicz (P-D) model adopting the same notation since exactly alike to the Biot one. 

The main assumptions of the theory are the same of the Biot theory. In particular: the 
solid phase of the system is perfectly elastic, homogeneous and isotropic; the liquid is a 
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compressible perfect fluid; the liquid viscosity is only introduced into a dissipation function 
to account for the sound wave attenuations due to the viscous friction between matrix and 
liquid; sound wavelengths are much larger than the heterogeneity dimensions; finally, all the 
transformations (displacements, strains, etc.) are considered infinitesimal in order to obtain 
a linear theory. 

In comparison with Biot model, two relaxation equations are added for the matrix and 
liquid temperatures and some thermo-mechanical coupling terms are inserted in all the 
equations to account for the coupling effects between the two temperatures and the two 
densities. 

Since in an isotropic medium TG experiment is sensible only to density and temperature 
changes, we focus our attention only to the motion equations of the matrix and liquid 
dilatations discarding the shear motions. 



Defining with u and U the displacements of solid and liquid phases J42| , the dilatations 
are defined by e = V ■ u and £^ = V • U. e and £ are related to the matrix and liquid density 
changes, 5 pi and 5p2, by the following identities: 

5pi = — Pie and Sp2 = — P2^ (7) 

The time evolution of the dilatations e and S and the temperature variations ST^ and 6Ti 
can be evaluated by the P-D model. In the appendix we report the detailed description of 
the P-D equations and all the expression of the coefficients. Few approximations suitable 
for the water- Vycor system can be made to this model (see appendix), in particular we can 
assume the local thermal equilibrium by which the liquid and matrix temperatures can be 
always considered locally equilibrated: i.e. T^ = Ti = T. After that, the linearized P-D 
equations in the Fourier q-space become: 

Piie + Pi2£ + Ke -£) + q^Pe + q^QS - q^RiST = 
Pue + P22S - b{e -£) + q^Qe + q^R£ - q^R25T = 

F6T + RiToe + R2To£ + q^k5T = (8) 

where the densities pij are related to the solid and liquid densities, ps and pi, by pu = 
(1 — ((>)fh — P12, P22 = (pPi ~ P12 and P12 = (1 — T)(f)pi where r is the tortuosity of the 
matrix 3|, |j]. The h{e — £^) is a viscous friction term arising from the relative motion 
between the two materials, being h = (p'^pi/kf), with pi the water dynamic viscosity and 
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ko the Vycor permeability. This term is, at the same time, the damping source for the 
acoustic waves and the hindrance to the hquid slow flow through pores. P, Q and R are 
the generalized isothermal elastic moduli, the coeflicients Ri are the gen eralized expansion 
coeflicient and the F coeflicient is a generalized speciflc heat [3|, \^ |43|, |4J] . k is the effective 
thermal conductivity. The complete deflnition of these coefficients is reported in appendix. 
These equations enable to calculate the g-component of the time dependent variables e, S 
and 6T, and thus, using the equations [H [5]and[7l to simulate the measured signal. 

The model of P-D model is an extension of Biot's model in the low frequency limit, i.e. is 
valid for motions at frequencies lower than the Biot characteristic frequency. This quantity 
is defined by /^ = fii/irpir^ with r the mean pore radius of the matrix. For the water- Vycor 
system the characteristic frequency is 26—80 GHz in the temperature range considered. This 
value has to be compared to the higher frequency experimentally excited which is almost 
0.6 GHz (the experimental frequency can be obtained by / = cq/27r where c = 4 Km/ s 
is the measured sound velocity and q = 1 fim~^ the exchanged wave vector). Thus, our 
experiment is testing dynamics at frequencies always lower then fc and the present model 
is theoretically suitable to describe our signals. Moreover, when the involved frequency 
is much lower than the characteristic frequency, the dependence on the tortuosity, r, of 
the solutions of Eq.s (lA14p . should be negligible. This condition has been verified in our 
case by the fitting procedure. As final consideration, we want to stress that the above 
equations, and consequently their solutions, do not depend explicitly on the pore dimension. 
Nevertheless, this parameter enters in the definition of the characteristic frequency, thus, 
fixing the applicability limit of the theory. 

B. TG response and signal 

The set of Eq.s (IHl) can be reduced to a first order differential equation system by intro- 
ducing two additional variables ^/'i = e and ?/'2 = ^- The set of equations ([8]) can, thus, be 
easily written in the compact form 

X(t) = -M ■ X(t) (9) 

where X = {e,£,ipi,ip2,^T) and M is a matrix of coefficients which can be easily obtained 
from Eq.s ([8]). Indeed we can solve the system by diagonalizing, with standard routines, the 
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non symmetric M matrix, i.e. M ■ V = V • D where D is diagonal and the resulting time 
solution is written as 

X(t) = V-exp(-D t)- [V-i-X(O)] (10) 

We note that it is possible to avoid the V matrix inversion, since Y = V^^ ■ X(0) is the 
solution of linear system V ■ Y = X(0). We want to stress that only the elements of the 
amplitude matrix, V, and the root matrix, D, need to be numerically calculated. Each 
element of the solution vector X(t) is a sum of two oscillating terms, cosine and sine rising 
from the two complex and conjugate roots of D, and three pure exponential terms deriving 
from the other three real roots of D. 

The temporal expressions of the dilatations and the temperature must be calculated 
considering that the electrostriction produces a non zero initial condition on ipi and 1^2 and 
the heating produces a no null initial condition on 6T: X(0) = (0, 0,'?/'i(0),'?/'2(0), 5T(0)). 
These solutions together with the expressions for the densities ([7]) and the expression for 
dielectric constant change ([5]) give the signal function: 

S"^{t) oc 6e{t) = Ae-'/^^cos{uJst) + Be-'/^^sin{uJst) + 

where the amplitudes A, B, C, D, and F, the acoustic parameters ts and us, and the 
time constants ri, T2, and T3 are functions depending on all the hydrodynamic parameters 
appearing in the starting equations and on the initial conditions and need to be numerically 
calculated at changing of the fitting parameters listed below. 

The oscillating terms in Eq. fllip accounts for the induced acoustic oscillations and the 
three exponentials for the remaining part of the signal. As previously discussed, the acoustic 
oscillations induced by the temperature grating are partially prevented. Anyhow, these ones 
have been equally considered in the fitting analysis and they proved to be very small at all 
the temperatures, around 5% of those induced by the electrostriction. 

Finally, Eq. (ITTl) with the addition of electronic response function and convoluted with 



the experiment instrumental function, yields the used fitting function 25 1. 
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VI. DATA ANALYSIS AND DISCUSSION 

The literature provides many of the water and Vycor data appearing in the P-D equations: 
the densities p^, pu the expansivities a^, a^ the water viscosity /i^, the specific heats Cp^, Cpi, 
the Vycor porosity 0, the transverse and longitudinal sound velocities of solid constituent 
for the calculus of Kg and the adiabatic sound velocity of water and the specific heat ratio 
7,, for the calculus of Ki. 

The other parameters were free fitting parameters. In table [T] we list the parameters 
locked to the literature values with the corresponding references. 

As reported in Sec. Illlt the solid constituent of Vycor 7930 is not pure silica but Vycor 



7913, a mixture of 96% of silica and 3% of boron oxide. Corning company |35| supplies us 
all the parameters of interest of Vycor 7913 only at room temperature. To overcome this 
problem, we made the hypothesis that the thermodynamic parameters of Vycor 7913 follow 



the temperature dependence as in silica j45| . Considering the restricted range of temperature 
analyzed, this is surely a valid approximation. 

For the calculation of the shear and bulk moduli of dry matrix, we take the matrix 
longitudinal sound velocity as free fitting parameter and the transverse sound velocity pro- 
portional to the longitudinal one as in the solid phase: Vt^m = 0.62VJ,m where the m and s 
indices refer to the matrix and solid materials and I and t to the longitudinal and transverse 



wave mode respectively 46|-|48 1 . 



As porosity, we use the value of 31%, which is the mean of values reported in the litera- 



ture 



12 



15 



49|. 



MM 

ce [7| ■w 



In reference [7| we proved that the Biot theory is not able to predict the acoustic damping 
times of TG data giving far overestimated values. The reason, probably, lies in the only 
damping mechanism considered in the theory, i.e. the energy dissipation due to the viscous 
friction between the matrix and liquid. Other damping effects, like the intrinsic dissipation of 
two materials are not included. P-D model does not introduce substantial modifications for 
the wave propagation and retains the same limitations of the Biot model. In order to avoid 
that the extraction of other interesting parameters could be affected by this mismatch, we 
multiply the oscillating terms of the solution (TTTj) by e"*/'^-'' to account for the real exponential 
decay of the acoustic waves. 

Finally, the free fitting parameters were: the three initial conditions '?/^i(0), 1^2(0) and 
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TABLE I: Water and Vycor parameters fixed in tlie fitting procedure to the values found in tlie 
corresponding references. 



Symbol 


Definition 


Values @ 20 °C 


Ref. 


Ps 


solid density 


2180 Kgm^"^ 


[35j 


Pi 


water density 


998 Kgm-"^ 


[50] 


as 


linear solid expansivity 


7.5 X 10^^ A'~i 


[35j 


Oil 


volume water expansivity 


2.07 X 10-^ K-^ 


[50] 


Pi 


water viscosity 


10-3 Pa ■ s 


[50] 


^ps 


solid specific heat 


739 JK-^Kg-^ 


m 


Cpi 


water specific heat 


4184 JK-^Kg-^ 


[50] 


4> 


porosity 


31% 


[12, 15, i6, 49] 


Vl,s 


solid longitudinal sound velocity 


5779 ms-i 


[35,45] 


Vt,s 


solid transverse sound velocity 


3580 ms"^ 


[35,45] 


Vi 


water adiabatic sound velocity 


1483 ms-^ 


[52] 


It 


water specific heat ratio 


1.0065 


[52,53] 



(5T(0), the relative amplitudes E and y4i2, the longitudinal sound velocity of the matrix 
Vi^mi the acoustic damping time ta-, the permeability ko and, lastly, the effective thermal 
conductivity k. This last parameter could to be fixed to the value given by the expression 
used in the the P-D model k = {1 — (j))ks + (pki, where kg and ki are the solid and liquid 
thermal conductivity, respectively. According to our data fitting, this value turned out not 
to be correct and we were forced to leave it free. 

The used fitting function is able to reproduce the data in all temperature range analyzed. 
In the left panel of Fig. |5] we report in semi-log scale the fit-data comparison of the signal 
at 20 °C together with the corresponding discrepancy. In the right panel of the same figure, 
we show in linear scale the same fit-data comparison, together with the three different 
contributions appearing in expression (|5]). These are the matrix density contribution (green 
line), the water density (blue line) and the temperature one (orange line) contributions. The 
sum of these three contributions gives the shown fitting curve. As clearly visible, the water 
density change is the main term and it is the only one yielding the bump in the signal. Thus, 
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FIG. 5: Comparison between data (black line) and fit (red line) for the temperatm'e 20 °C with 
the relative discrepancy (left panel). In the right panel the three different contributions whose sum 
gives the fit of the datum at 20 °C are shown. These are the matrix density contribution (green 
line), the water density (blue line) and the temperature (orange line) contributions. The inset 
shows the same contributions in x-log scale. 

the long time dynamics is mainly related to the water density changes and this is clearly due 
to the great difference between the water and Vycor expansivities and to the stiffness of the 
matrix. The signal contribution induced directly by the temperature (third term of Eq. HJ 
(de/dTi) 6Ti(t)) is practically described by a single exponential decay, it does not show any 
acoustic oscillation and has an intensity weakly dependent on the temperature. We want 
to stress that, as regards the acoustic oscillations, the liquid and matrix move in phase (see 
the inset of Fig. [S]). Contrary, in the major part of the signal the water and matrix densities 
have opposite derivatives, i.e. when the liquid tends to expand the matrix tends to compress 
and vice versa. This implies an out-phase motion of the liquid and matrix as expected for 
the slow sound predicted by the Biot theory. 

Once we extracted all the parameters for a given temperature, we can use them to simulate 
the time dependence of an isolate physical observable. It is interesting to study the temporal 
behavior of the liquid pressure inside the pores. This can be obtained from the constitutive 
equation for the pore fluid [28|: 



Spit) = --[Qeit) + R£{t) - (R, + R2)STit)] (12) 

Figure [6] shows the temporal evolution of the pore pressure change calculated with the 
parameters of the fit at 20 °C. It is interesting to note that the liquid pressure change is 
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FIG. 6: Temporal evolution of the pore pressure obtained from the constitutive equations for the 
pore fluid (Eq. [T2]) using the parameters extracted by the fit of the datum at 20 °C. The pressure 
shows two relaxation modes, a fast one which drives the acoustic oscillations (inset in the figure) 
and a slow one which drives the liquid motion inside the pores. 

never equal to zero (apart when it changes sign) during the signal evolution. This behavior is 
different from what appears in a bulk liquid sample where the pressure equilibrates after the 
vanishing of the acoustic oscillations. In the case of confined water, the liquid expansion due 
to the heating is partially prevented by the stiffness of the matrix which exerts a pressure on 
the liquid. At first, the liquid will change its density via the outflow from the heated pores 
to nullify this pressure. At later times, owing to the vanishing of the thermal grating for 
thermal diffusion, the matrix will exerts a negative pressure on the liquid which backflows 
to the pores to equilibrate again the pressure. Thus, density changes related to the liquid 
flow take place on time scales much longer and the liquid pressure is different from zero at 
all the times. Finally, the liquid pressure shows two relaxations, a fast one which drives 
the acoustic oscillations (the fast Biot wave) and a slow one which drives the liquid motion 
inside the pores (the slow Biot diffusive wave). 

VII. FITTING RESULTS 

In this last section, we shall show and discuss the main flndings derived by the flt- 
data analysis. We will start with the results of acoustic wave propagation after that we 
will analyze the features connected with liquid viscous flow and flnally with the thermal 
diffusion. 
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FIG. 7: Longitudinal sound velocity of vycor as a function of temperature. 

As already said in previous section, many of the parameters appearing in the equations of 
P-D model have been kept fixed to the literature values. These values are, clearly, measured 
for the bulk materials. For the case of water the assumption to consider the thermodynamic 
and dynamic parameters of confined water equal to the bulk ones is not so obvious. The 
strong hydrogen bonds of water molecules with the silanol groups of the inner pore surfaces 
of Vycor could affect parameters like the density, specific heat, viscosity, thermal expansivity, 
etc. For example, it has been found by molecular simulation 5J] that the density of water 
confined in Vycor is around 11% lower than the bulk value. Yet, the specific heat has 
been measured 1 — 2% higher than in the bulk 55|. Contrary, measurements of viscosity 
of water confined in nanometer films proved that this parameter remains close to the bulk 



value 



56 



571]. Nevertheless, the data about the thermodynamic parameters of water in 



confined state are few and, in particular, no data about their temperature evolution exist. 
For this reason, we have chosen to fix the values of these parameters to the bulk ones as it 
has been done so far in all the experimental study on confined water. 

In Fig. [7] we report the temperature behavior of the longitudinal sound velocity of Vycor 
yi,m- We were forced to leave this parameter as free fitting parameter in order to have a 
good fit of acoustic oscillations at all the temperatures. These values together with the 
transverse velocity values, Vt^m = 0.62Vi^rn, and the values of matrix density, give the matrix 



stiffness values which differs from the those measured by static experiments [15|, |16|, |35| . But 
our longitudinal sound velocities are in reasonable agreement with the values measured by 
Brillouin scattering experiments by Levelut and Pelous 48 1. 
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FIG. 8: Measured acoustic attenuation times (circles) compared with the values predicted by the 
model of Pecker and Deresiewicz (continuous line). 

In Fig. [8] we show the measured acoustic attenuation times compared to the prediction of 
the P-D model. We see that the times predicted by the theoretical model are much longer 
than the experimental ones and moreover their temperature behavior show an opposite 
trend. As stated before, P-D model does not introduce substantial modifications for the 
wave propagation and retains the same limitations of the Biot model. As just reported 
in [TJ , the only damping mechanism considered in the theory (the energy dissipation due to 
the viscous friction between the matrix and liquid) is not surely the main damping source in 
Vycor filled with liquids at our TG frequencies. Other damping effects should be included, 
beginning with the dissipation phenomena intrinsic of the Vycor and water. Considering 
the very low sound attenuation in bulk water, the damping could be mainly due to the 
sound damping in Vycor. Indeed, the underestimations of the damping mechanism has 



been already reported in the literature 581-160 1. 



In Fig. |9] we show the measured permeability ko as a function of temperature, ko does 
not show, within the experimental errors, a temperature dependence and has a mean value 
of ~ 0.126 nrn^. The permeability is related to the viscous Darcy coefficient b = (p'^fii/k^ 
and to the hydrauhc diffusivity coefficient p]: Djj = {PR — Q'^)/[b{P + R + 2(5)], being 
P, Q and R the generalized elastic moduli defined in the appendix. In fact, the P-D model 
is an hydrodynamical theory that supposes that the flow of liquid through the pores is a 
Darcy flow (i.e. that the flow is laminar and the fluid velocity at pore wall is zero). This is 
assumed in the deflnition of the dissipative force b{e — S) and, in particular, in the expression 
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FIG. 9: Permeability of water in Vycor extracted from the fits. 

of b. In this model, the permeability of the porous material, fc^, is a quantity depending 
only on the geometrical characteristics of the porous material, porosity, mean pore size 
and tortuosity. Indeed the kn values extracted by the fit does not show any appreciable 
temperature dependence. 



The value of the measured permeability can be compared with that of reference 15|, 
161 ] [ki) = 0.065 nm^) and with that calculated following the expression of the Poiseuille 
permeability for a porous medium 6l|: k^ = (l/8)r^0/r where we recall r and r are the 



mean pore radius and the tortuosity of the medium. By inserting the value of 31% for 



the porosity and of 2 nm 
a tortuosity value of 2-4 



'or t 



le mean pore radius, supplied by the Corning company, and 
13l . Iisl . [l6|, we find a permeability of 0.04 — 0.08 nrv?. Thus, 



our kr) fitting values do not agree with the previous reported measurements. This could 
be ascribed to a mean pore size larger than what stated by the company. Frequently in 
literature, measured values of porosity and mean pore size have been found different from 

nn 

the usual Corning values. For example, in ref. |15l. 1161]. it is reported a mean pore radius of 
2.68 nm in [ij] of 2.3 nm. Our results would be compatible with a mean pore radius around 
3 nm. 

In Fig. [in] we report the temperature behavior of the measured effective thermal conduc- 
tivity, k. We report also the values of the thermal conductivities calculated according to 
the weighted on porosity arithmetic mean (red line), k = {1 — (f))ks + (ph, and harmonic 
mean (blue line), 1/k = {1 — (f))/ks + (p/h, being kg the solid thermal conductivity equal 
to 1.29 Wm'-'^K-'^ @ 20 °C and ki water thermal conductivity equal to 0.60 Wm'^K'^ @ 
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FIG. 10: Effective thermal conductivity as a function of temperature compared with the maximum 
theoretical values (red line) and the minimum ones (blue line). 



20 °C SOj. These expressions of the effective conductivity refers to a medium composed by 
parallel layers of solid and water where the heat transfer goes along the layer direction (par- 
allel conduction) or goes perpendicular to the layer direction (series conduction) 20|. The 
value of these two values can be considered as the two limiting values of the k parameter, 
see also appendix. The values obtained by our fitting is in substantial agreement with the 
harmonic mean. 

Finally, we want to show the results obtained for the amplitude parameters E and A12. 
These are shown in Fig. [TTJ Both the parameters do not show, within the experimental error 
bars, a sensible temperature dependence. In the same figure, the E values are compared 
with the ones obtained for the bulk water |2J]. Unfortunately, we have not enough sensibility 
to estimate if and how the confinement affects the photothermal effect of water. Anyway, 
we want to underline that both the matrix and temperature contributions are absolutely 
necessary to describe correctly the data. 

Considering the complexity of the fitting model, the agreement of the P-D parameters 
with the literature data is reasonable even if it is not complete. Correction on some fixed 
parameter could lead probably to a better agreement. In particular, possible changes of the 
water parameters due to the confinement should be taken in to account. Nevertheless, the 
fact that the model is able to reproduce the data at all the temperatures with a substantially 
small number of free parameters catching the main temperature dependencies and the order 
of magnitude of the coefficient values, proves the overall validity of the hydrodynamic laws 
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FIG. 11: Temperature behavior of the parameters E (blue chxles) and A12 (red squares). Green 



triangles refer to the E values measured in the bulk water 



2M- 



in describing an nano-heterogenous system. 



VIII. CONCLUSIONS 



In summary, heterodyne detected transient grating experiment has been applied to study 
the relaxation dynamic of water confined in Vycor 7930. We acquired HD-TG data in the 
temperature range 20 — 90 °C at the q- vector of 1 fim^^ and for the range of g = 0.63 — 2.5 
^m~^ at 30 °C . The HD-TG experiment enables to investigate in a single data the damped 
acoustic oscillations taking place at short times, the mass/density dynamics and the thermal 
diffusion. The HD-TG data has been analyzed using the hydro dynamic model of thermo- 
poroelasticity on porous media of Pecker and Deresiewicz J28[. This hydrodynamic model 
gives a very interesting insight of the experimental results enabling a consistent charac- 
terization of the transport processes as they appear in the experimental signal. We need 
to redefined few parameters of the model in order to reproduce correctly our experimental 
data. In particular the acoustic damping rate that the P-D model, as well as the Biot theory, 
underestimate drastically and the effective thermal conductivity that can not described as 
a simple weighted mean. The simulation of HD-TG data according to the P-D model shows 
that acoustic propagation, taking place on the fast time scale, is due to an in-phase motion 
of the solid and liquid components of the system. The mass/density transport is due to the 
flow of the confined liquid that results to be coupled to the heat diffusion. A peculiar flow 
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and back-flow of water inside tlie nano-pores forced by thermal processes has been reported. 
During such transport processes the solid and the liquid are moving with opposite phases. 
We fitted our data using the P-D model, fixing as much as possible the water- Vycor param- 
eters to the known literature values. The values of the free parameters have been extracted 
from our data using a best fit procedure. The P-D model enables a valid fit of our data 
for the whole investigated time windows, even if there is not a complete agreement of some 
fitting values with the values reported in the literature. In our opinion this disagreement is 
not to be ascribed to a fail of the hydrodynamic laws but to an oversimplified definition of 
the solid and liquid basic parameters. 

We would like to stress that HD-TG experiments turns to be able to measure the complex 
dynamic processes, taking place in a nano-heterogenous system, over a very broad time 
windows and a relatively simple hydrodynamic model seems appropriate to describe all 
the transport and dynamic phenomena measured. Our study gives a further experimental 
confirmation of the fact that water flow in Vycor is well described by hydrodynamic laws in 
spite of the nanometric dimension of pores. 
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Appendix A: The Pecker and Deresiewicz hydrodynamic model 

According to the Pecker and Deresiewicz model the linearized equations describing the 
time evolution for the dilatations e and £ and the temperature variations 5Tm and 5Ti 



are |28|: 



Pu'e + Pi2£ + h{e -8)- PV'e - QV^E + RuV^STm + Ri2^^5Ti = (Al) 

Pue + P22£ - h{e -£)- QV'e - RV^£ + i?2i V'^T^ + i?22 V'cJT; = (A2) 
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FnStrn + Fi25Ti + K{5Trn - m) - (1 - 0)fc.V25T^ + R^{r^e + i?2iTo^ = (A3) 



F2i(5t^ + F^^m - K{6T„, - 6Ti) - <PhVHTi + R^^T^e + R22To£ = 



(A4) 



in which 

- the densities pij are related to the sohd and hquid densities, ps and p;, by pn = 
(1 — 0)/?9 — pi2, P22 = 0Pi — P12 and P12 = (1 — t)(J)Pi where r is the tortuosity of the 



matrix 



U- 



h = (fP'pi/ko-, with /i; the water dynamic viscosity and ko the Vycor permeabihty. The 
h{e — £^) is a viscous friction term arising from the relative motion between the two 
materials. This term is, at the same time, the damping source for the acoustic waves 
and the hindrance to the liquid slow flow through pores. 

The P, Q and R coefficients are the generalized isothermal elastic moduli defined in 
the Biot theory [sl, U], they can be related to the isothermal bulk modulus of liquid 
Ki, the bulk modulus of solid Ks-, the bulk modulus of matrix Krn and to A^ which is 
the isothermal shear modulus of matrix. In particular we have 43|, |44| 



P 



Krn 
Ks 



Krn -LAKs 

Ks ^"f" Kl 



+ -N 

3 



Q 



R 



1- 



^)4>Ks 



1- 



Krn _L(/,i^ 
Ks ^^Ki 



Kn 
Ks 



l-S- ^ + (h^ 



Ki 



(A5) 

(A6) 
(A7) 



with the matrix porosity. 

We recall that A';, A'^, K^ and A^ can be expressed as functions of the density and 
the longitudinal and transverse sound velocities 



Ki = p.Vlhr , K, = psiV^l - 4/3VI) 



K„ 



Pm(V^' -4/3V;' ), iV = P,n^M 
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The coefficients Rij are the generahzed expansion coefficient defined by following equa- 



tions 



43 



4^: 



Ru = a,{3P - AN) + aisQ 
Ri2 = asi{3P - 4N) + aiQ 

R21 = 'iasQ + aisR 

R22 = SasiQ + aiR (A8) 

where as is the linear expansion coefficient of solid and ai is the volumetric thermal 
expansion coefficient liquid, asi and ais are the thermo-elastic coupling coefficients. 

The Fij coefficients the are generaliz ed spec ific heats at constant volume and they are 



defined by the following expressions 43|, |4J] : 



Fii = (1 - (t>)psCps - To(3a,i?ii + aisR2i) 
F12 = -To(3asi?i2 + aisR22) 
F21 = -To{3asiRii + aiR2i) 
F22 = (ppiCpi - To{3asiRi2 + OiiR22) (A9) 

where Cps and Cpi are the isobaric specific heats of solid and liquid and Tq the equi- 
librium temperature of the whole system. 

- K is the coefficient of interphase heat transfer. 

- kg and ki are the thermal conductivities of solid and liquid. 

Now some approximations can be made to simplify the equations. As already discussed 
in section IIV| we can assume the local thermal equilibrium between the two phases. This 
implies the assumption of a if coefficient very high in the equations and thus T^ = Ti = T. 
An other hypothesis is that the thermo-elastic coupling coefficients, agi and ais can be 
considered negligible. This is a plausible approximation considering the very low thermal 
expansion coefficient of Vycor and its stiffness. In fact, agi is the strain in the matrix due to 
an unit change of temperature in the liquid and ais is the dilatation of the liquid due to an 
unit change of temperature in the matrix. We expect that the former is negligible due to the 
stiffness of the matrix {Km = 13.4 GPa > Ki = 2.2 GPa @20 °C) and the latter is negligible 
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due to the low expansivity of the matrix (Sa^ = 2.25 x lO"*" K~^ « ai = 2.07 x 10""^ K~^ 
@20 °C). After these approximations and summing together Eq.s (1A3P and f lA4p in order 
to obtain a single equation for the temperature T, the system of equations (lAip becomes 

Piie + puS + b{e - £) - PV\ - QV^S + RiV^6T = (AlO) 

Pue + p22^ - b{e -£)- QV^e - RV^£ + R2V'^6T = (All) 

F6T + RiToe + i^aTo^ - W^6T = (A12) 

where 

Ri = as{3P - AN) + aiQ 

F = ^ Fjj = (1 - 0)psCp, + (ppiCpi - To(3asi?i + a^i^a) 

k = {l- (f))ks + (j)ki (A13) 

P-D model, in its original form, predicts an effective thermal conductivity, k, defined by 
an arithmetic means of values of the two constituent phases, weighted on porosity: k = 
(1 — (f))ks + (f)ki. This kind of mean is used also in the definition of the effective volumetric 
heat capacity appearing in the expression of F coefficient. Since the effective specific heat is 
independent from the medium morphology, this definition appears appropriate and correct. 
Whereas, it is well known that the thermal conductivity is generally dependent by the 



medium morphology J20|], in the present sample by the geometric characteristics of the 
porous material. The expression reported in eq. IA13I refers to a medium composed by 
parallel layers of solid and water where the heat transfer goes along the layer direction 
(parallel conduction). Clearly, this picture represents a very particular case which, anyway, 
defines a maximum limit for the effective conductivity. The minimum limit is, instead, 
obtained with the harmonic mean, 1/k = (1 — cfy/kg + 0/^;, which describes a composite 
medium of parallel layers where the heat transfer goes perpendicular to the layer direction 
(series conduction). In a generic porous material, where the tubules are interconnected and 
randomly oriented in all directions, the effective thermal conductivity should assume values 



within these two limits 



20 



62 



63| 
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Since the TG signal is directly related to a g-component of spatial Fourier transform of 



the dielectric constant change 25|, the next step is writing in the g-space the system of 
equations (lA12p : 

Pii'e + Pi2£ + h{e -£) + q^Pe + q^Q£ - q^RiST = 
Pi2e + P22S - Ke -£) + q^Qe + q^R£ - q^RiST = 

F6T + RiToe + i^sTo^ + q^kST = (A14) 

where, now, e, S and 6T refer to a g-component of the their spatial fourier transforms. 
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